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We analyze the energetic and dynamical properties of bright-bright (BB) soliton pairs in a binary 
mixture of Bose-Einstein condensates subjected to the action of a combined optical lattice, acting as 
an external potential for the first species, while modulating the intraspecies coupling constant of the 
second. In particular, we use a variational approach and direct numerical integrations to investigate 
the existence and stability of BB solitons in which the two species are either spatially separated 
(split soliton) or located at the same optical lattice site (overlapped soliton). The dependence 
of these solitons on the interspecies interaction parameter is explicitly investigated. For repulsive 
interspecies interaction we show the existence of a series of critical values at which transitions from 
an initially overlapped soliton to split solitons occur. For attractive interspecies interaction only 
single direct transitions from split to overlapped BB solitons are found. The possibility to use split 
solitons for indirect measurements of scattering lengths is also suggested. 

PACS numbers: 67.85.Hj, 03.75.Lm, 03.75.Kk, 67.85.Jk 


I. INTRODUCTION 

Bose-Einstein condensates (BECs) are fascinating 
tools for simulating different physical systems. Advanced 
laser technology and its successful applications to ultra¬ 
cold atoms have enabled us to engineer potentials of dif¬ 
ferent geometries. A well-established technique consists 
in creating a linear optical lattice (LOL) by interfering 
pairs of counter propagating laser beams [1]. On the 
other hand, laser beams can also be used to vary atomic 
interaction periodically in space with the help of optical 
Feshbach resonances [2]. Periodically modulated atomic 
interaction leads to a nonlinear optical lattice (NOL). 
The LOL has been used to investigate different physi¬ 
cal phenomena in condensed matter physics, including 
Bloch oscillations laii], generation of coherent atomic 
pulses (atom laser) [^, dynamical localization la m, 
Landau-Zener tunneling [HHin] and superfluid-Mott tran¬ 
sitions El- 

Interatomic interaction in BECs gives rise to a non¬ 
linearity which permits localized bound states to remain 
stable for a long time, due to the balance between the 
effects of nonlinearity and dispersion. In the presence 
of a LOL, the interplay between lattice periodicity and 
interatomic interaction was shown to induce modulation 
instabilities of Bloch wavefunctions near the edge bands 
[12], leading to the formation of localized excitations 
with chemical potentials inside band gaps, the so-called 
gap solitons (GSs). These excitations have been investi¬ 
gated both for continuous BECs, in one-dimensional [ni¬ 
di] and multi-dimensional [T8]-[20j settings, and for BEC 
arrays in the presence of attractive and repul¬ 

sive interactions. NOL can also support special kinds of 
solitons both in ID [n] and in multi-dimensional set¬ 
tings in combination with LOL [24l [25]. NOLs have 
been used to avoid dynamical instabilities of gap-solitons 


and to induce long-lived Bloch oscillations [26], Rabi 
oscillations m and dynamical localization |28j in the 
nonlinear regime. Eor comprehensive reviews on single¬ 
component BECs in linear and/or nonlinear optical lat¬ 
tices see [531 153H3T] . 

On the other hand, the analysis of the physical proper¬ 
ties of binary mixtures of condensates still displays open 
issues, and represents an interesting research topic [31]- 
138] . In the past years some work has been done on the 
stability and dynamics of binary BEC mixtures with both 
components loaded in LOLs [33] or in NOLs [JD] or com¬ 
binations thereof [311144] . 

However, BEC mixtures with one component loaded 
in a LOL and the other loaded in a NOL have not been 
investigated, to the best of our knowledge. This setting 
is particularly interesting because it may support new 
types of matter waves, due to the interplay between the 
different types of OL and the intrinsic nonlinearities. In 
particular, in absence of any interaction (e.g. with all 
scattering lengths tuned to zero), the spectrum of the 
component in the LOL displays a band structure, while 
that of the other component has free-particle features. 
It is known that for attractive intraspecies interactions, 
uncoupled mixtures will feature localized states. In this 
situation one can expect that a rich variety of bound 
states can be formed once the interspecies interaction is 
switched on. 

The aim of the present paper is to study localized mat¬ 
ter waves of binary BEC mixtures with one component 
loaded in a LOL and the other in a NOL. In particular, we 
concentrate on localized states which have chemical po¬ 
tentials of both components in the lower semiinfinite part 
of the spectrum. We call these states bright-bright (BB) 
solitons, or also “fundamental” solitons, because when in¬ 
traspecies scattering lengths are both negatives (the case 
investigated in this paper) they coincide with the ground 
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state of the system. We show that BB solitons can be 
classified according to the distance between the lattice 
sites where centers of their components densities are lo¬ 
cated. Denoting these distances by nL, n = 0,1,2,..., 
with L the spatial period of the lattices (assumed to be 
the same for both LOL and NOL), the n = 0 and n 0 
families are referred to as overlapped and split BB soli¬ 
tons, respectively. The existence and stability of these 
solitons are investigated both by a variational approach 
(VA) for the mean-field two-component Gross-Pitaevskii 
equation (GPE), and by direct numerical integrations of 
the system. In particular, the dependence of the exis¬ 
tence ranges of BB soliton pairs on the interspecies inter¬ 
action parameter, 712 , is investigated. As an interesting 
result, we find that one can pass from one soliton family 
to another by simply changing the strength of the inter¬ 
species interaction. In particular, starting from an over¬ 
lapped (n = 0) BB soliton one finds a series of repulsive 
critical values of 712 at which the transition from the n- 
to the n + 1-split BB soliton occurs as 712 is adiabatically 
increased away from the uncoupling limit (712 = 0). On 
the contrary, for attractive interspecies interaction only 
direct transition from split to overlapped BB solitons are 
possible. Since critical values at which transitions oc¬ 
cur depend on physical parameters of the mixture, these 
phenomena suggest that split BB solitons could be used 
for indirect measurements of scattering lengths in real 
experiments. 

The paper is organized as follows. In Section II, we in¬ 
troduce the mean field equations for the coupled system 
and envisage a variational study for stationary localized 
states. We examine the linear stability of these states 
for attractive and repulsive intercomponent interaction. 
In Section III, we introduce a time-dependent variational 
approach, with Gaussian trial solutions, to study differ¬ 
ent classes of BB soliton pairs. The stability of split and 
overlapped families of soliton pairs is checked by numeri¬ 
cal integration of the mean-field equations. In Section IV, 
a numerical routine is employed to understand the role 
of interspecies interaction in the splitting mechanism, for 
both attraction and repulsion between different species. 
Finally, in Section V we make concluding remarks. 


II. ANALYTICAL FORMULATION 

Throughout this paper, we shall consider a quasi-one- 
dimensional binary mixture of BEGs, in which the trans¬ 
verse motion is frozen into the ground state of a tight 
transverse trapping potential, with trapping frequency 
a;_L. The mean-field dynamics of a mixture in which the 
two species’ particles have equal mass m is modeled by 
the coupled GPEs [JS] 

-|-2ai2^W_L|d>3_jp J dtj, (1) 


where j = 1,2 is the species index, Vj’s are the exter¬ 
nal trapping potentials, Uj’s the intraspecies scattering 
lengths (which generally depend on position) and 012 
the interspecies scattering length. The wave functions 
are normalized to the numbers of particles 

A/; =|ds|vE-,f. (2) 


Since our system is subject to an external potential 
proportional to cos{2kLs) generated by two counterprop- 
agating laser beams, the inverse wavenumber and the 
recoil energy Er = (Iifci)^/ 2 m provide natural units for 
length, energy and time |29] . To simplify the notation, 
we introduce the adimensional quantities 


X 

= kLS, 


V, 

Vi 



2Er’ 

N, 

huj\ 

— 


2Er 


2Er 


hu!± 




2ErkL 

Afj, ^j:=2ajkL, 712 := 2012 ^^, (3) 


yielding the GPEs 

* ^ ^ ^ ^ I I ^ I I ’ 

(4) 

and the constraint 



(5) 


(Note how the variables Nj , which will be used through¬ 
out this article, coincide with the actual numbers of par¬ 
ticles Afj only up to a factor.) In the physical case of 
interest, only the first species is subject to an external 
lattice potential: 


Vi{x) = Vbi cos(2a;), V 2 {x) = 0. ( 6 ) 

The loose longitudinal harmonic trapping will be ne¬ 
glected, since we will focus on states that are localized 
over a few lattice sites. As for the interspecies coupling 
constants, we shall assume that the interaction of the 
second-species particles depends on the lattice modula¬ 
tion: 


7 i(a:) = 7 i, 72 ( 0 :) = 72 -f Vb 2 cos(2a:). (7) 

To allow the existence of bright solitons, which are sta¬ 
tionary states with both species localized, the average 
coupling constants are assumed to be negative (so that 
I 72 I > | 4 o 2 |)- Matter-wave bright solitons have also been 
observed experimentally in trapped systems m- Since 
the presence of the linear and nonlinear lattice poten¬ 
tials in @ is an obstruction to finding exact bright soli¬ 
ton solutions, our study will be based on a reasonable 
variational approach, with a subsequent numerical test. 
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FIG. 1 . Left panel: dependence of the minimum energy of the 
system, within the Ansatz vs the number of atoms A^i, 
with Ni + N2 = 3 . Right panel: difference between chemical 
potentials = /ri — /22 vs Ni. The parameters are fixed to 
71 = 72 = —l,Vbi = — 0.5 and V02 = — 0 . 25 . Energies and 
chemical potentials are in units of 2 Er. 


A. Stationary solutions: overlapped solitons 




FIG. 2 . Existence curves in the fj,i — N plane {N = Ni = N2) 
for BB solitons of the GPEs with attractive (712 = — 0 . 5 , 
left) and repulsive (712 = 0.1, right) interspecies interac¬ 
tion. Other parameters are fixed as 71 = 72 = —l,Voi = 
— 0 . 5 , V02 = — 0 . 25 . Square and triangle lines denote numer¬ 
ical results obtained from Eqs. by imaginary time relax¬ 
ation method and refer to the hrst and the second compo¬ 
nent, respectively. Dashed and dot-dashed lines represent the 
corresponding results obtained from the minimization of the 
energy functional. 


We are interested in the stationary solutions of the cou¬ 
pled GPEs Q. The form yields 

the stationary GPEs 

/ 1 9^ \ 

(,“ 2 ^ + +7i2|<^3-jf = 0, 

( 8 ) 

where the external potentials and coupling constants are 
given, respectively, by (|^ and Q. Each stationary 
state is characterized by the chemical potentials (/ii, /i 2 ), 
which are fixed by the normalization conditions. 

Let us assume that the intraspecies interactions are 
attractive ( 7 ^- < 0). Moreover, we shall focus on the case 
^ 01)^02 < 0 : in this situation, due to the (linear and 
nonlinear) trapping mechanisms, density profiles peaked 
around the points where cos( 2 a:) = 1 are energetically 
favorable for both species. In order to investigate the 
features of BB soliton pairs, we choose a Gaussian trial 
solution 


function of the soliton width: 


E = dx 


sEI 

. i=i,2 




dx 




+ 7 i 2 |</'in</> 2 p + Vi{x)\cj)i 


1 / V^-^1 \/^N2 /T^rr 

_l_ V02N2 ^ 71-^1 72-^2 I12N1N2 \ 

0-2 0,1 O2 ^a\+ 0,2 ) 

( 11 ) 

The optimal width values (di(Ai, A^2), d2(iVi, A2)) are 
determined by 


dE 

doj 


= 0 for j = 1 , 2 , 


( 12 ) 


4>j{x) = Aj exp [—a:^/2a|] . (9) 

Since the amplitudes Aj and the widths Oj are bound by 
the normalization conditions 

Nj = J dx\c/)j{x)\^ = \^ajAj, ( 10 ) 

the functions ([^ have only one free parameter. More¬ 
over, this class of trial solutions fits overlapped BB soli¬ 
tons, with the peak of their densities sitting at the same 
position, say, at a: = 0. Since, due to attractive inter¬ 
species interactions, the superposition of densities lowers 
the energy of the system, we expect the most energeti¬ 
cally favorable soliton pair to be overlapped. 

At fixed numbers of particles, the Gross-Pitaevskii en¬ 
ergy functional for (j)j in the class can be viewed as a 


and fix the trial ground state. The corresponding energy 
will be denoted by 


Erain{Ni,N2) — E\a-=aj(Ni,N2)- 


(13) 


The chemical potentials p.j = dE^in/dNj can be used 
to test the linear stability of the ground state solution 
through the Vakhitov-Kolokolov criterion [46]. 

Relevant properties of the overlapped solitons can be 


inferred from the energy functional in Eq. (11) and the 


chemical potential. If one keeps constant the total num¬ 
ber of atoms N = Ni + N 2 , the change in the energy 
Eraini^i, N 2 ) of the trial ground state with Ni (or equiv¬ 
alently N 2 ) can be analyzed. Let us fix for definiteness 
7 i = 72 = -1, 712 = -0.5, Vbi = -0.5, V 02 = -0.25. 
Throughout the paper, numbers Nj of order one will be 
extensively used: recall that they are related to the actual 
numbers of atoms Nj by the factor 2Er/hujj_ thru Eq. (|^. 
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FIG. 3 . Profiles of BB solitons with Ni = N 2 = N corre¬ 
sponding to the A'’ = 1 points on the cnrves of Fig. for 
attractive (712 = — 0 . 5 , left panel) and repnlsive (712 = 0 . 1 , 
right panel) interspecies interaction. Blue (dashed and dot- 
dashed) and red (continuons and dotted) curves refer to first 
and second component, respectively. In both panels continu¬ 
ous and dashed lines refer to numerical solutions of Eqs. Q 
found by relaxation method, while dotted and dot-dashed 
lines to variational analysis. 


In an experiment with ^Li atoms (m ~ 1.2 x 10 kg), 
with a transverse trapping potential a;_L ~ 27r x 700 Hz 
EZlllH] and a laser wave number fci = 5 x 10 ^ m ^, the 
conversion factor is of the order of lO"*. 

In Fig. we set the (scaled) number of particles to 
N = 3, and study the behavior of the minimal energy 
and the difference in chemical potentials of the soliton 
pair. The evident asymmetry in the plot of the energy 
E{Ni,N — Ni) (left panel) with respect to = N/2 is 
related to the inhomogeneity of the lattice potentials and 
the self-interactions for the two components. Moreover, 
the difference in chemical potentials A/i = /ii —/r 2 (right 
panel) can lead to changes in the numbers of particles if 
the system is in contact with a particle reservoir or if a 
transition mechanism between the two species is present. 

Existence curves in the ^—N plane are plotted in Fig.[^ 
for the case of equal attractive intraspecies interactions 
and for both attractive (left panel) and repulsive (right 
panel) interspecies interaction. To reduce the number of 
parameters we have considered the case Ni = N 2 = N. 
The dotted lines refer to the VA results obtained from 
the numerical minimization of the energy ( [TT| ), while the 
solid lines represent the corresponding curves obtained 
from numerical relaxation method and self-consistent di- 
agonalization of the stationary GPEs in Eq. (HI dsmo]. 
From this plots one can see that d^j/dN < 0, both for 
the VA and for the GPE curves, which implies, accord¬ 
ing to the Vakhitov-Kolokolov criterion, that BB solitons 
are linearly stable. This result is also confirmed by di¬ 
rect numerical time integrations of the two component 
GPEs 0 (not shown here). 

From Fig. one can see that, while in the attractive 
case the agreement is quite good for a wide range of N 
(this is true also for relatively large values of 712 ), in the 
repulsive case deviations of the VA and numerical curves 
are larger. This discrepancy is due to the fact that the 
Gaussian Ansatz becomes less accurate in the repulsive 


case, in which the interspecies interaction reduces the sta¬ 
bility of the overlapped configuration, as one can see from 
Fig. where VA and GPE solitons profiles are compared. 
However, the accuracy of the VA result increases as posi¬ 
tive 712 values are decreased towards the uncoupling limit 
712 = 0. Also note the existence of points where the 
chemical potential curves intersect, both for attractive 
and repulsive interspecies interactions. At these points, 
BB solitons, having the same number of atoms and the 
same chemical potentials (related to their width), will 
have equal VA profiles for the two components. Despite 
the discrepancy in the location of the intersection point, 
equality of profiles is well confirmed by numerical GPE 
results. 


III. SPLIT BB SOLITONS AND DYNAMICAL 
PROPERTIES 


In the previous section, the choice of Gaussian trial 
wave functions ([^ aimed at studying overlapped BB soli¬ 
ton configurations, with the peaks of the two density 
profiles coinciding at a: = 0. Due to attractive inter¬ 
species interaction, this configuration is expected to be 
the lowest-energy BB soliton pair. It is possible however 
to extend the analysis to split BB solitons, in which the 
centers of mass of the two species do not coincide. 

Let us initially consider the uncoupled limit 712 = 0. 
In the case Vbi < 0 and Vb 2 < 0, one expects an infinite 
set of degenerate energy-minimizing BB solitons, since 
the centers of mass of each species can be located at any 
point Xoj such that cos( 2 a;oj) = 0 , regardless of the other 
species’ density profile. These minimizing configurations 
can be classified in families 


BBn(^i 2 = 0), with n e N, 


(14) 


according to the absolute distance between the centers of 
mass: 


Ai := 1x02 - xoil G 


1 


1 


"■2r’r+2 


(15) 


Clearly, the energy of the BB soliton configurations is the 
same for all families, En{'yi 2 = 0) = Eq. When the effect 
of the interspecies coupling 712 can be treated as a small 
perturbation, one expects the existence of stationary 
Gross-Pitaevskii solutions close to the ones at 712 = 0, 
which can still be classified in families BBn{ji 2 ) accord¬ 
ing to the criterion (15). However, if 712 < 0, the in¬ 


terspecies attraction will break the energetic degeneracy 
in favor of the overlapped configuration BBq. Thus, the 
split solitons in BBn>o{ji 2 ) become metastable. While 
configurations with a very large distance between the two 
species are almost unaffected by interspecies interactions, 
larger values of I 712 I weaken the (local) stability of soli¬ 
tons with small n, since attractive interactions can give 
a sufficient amount of energy to overcome the (linear or 
nonlinear) potential barrier and reduce the distance be¬ 
tween the centers of mass. Thus, one expects a critical 
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- 0.06 - 0.04 - 0.02 
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0.9293 0.9295 

{aun) 



2.32 - 
2.3 ■ 

L 2.28 ■ 

2.25 - 
2.23 


712 


0.9291 0.9294 0.9297 

(ai,n) 


712 

n 

Ax 

ffli 

02 

nmin 

0 

0 , 1,2 

0 , TT, 27r 

0.929 

2.229 

-0.177 

-0.25 

0 

0 

0.910 

1.472 

-0.249 

-0.25 

1 

no local minimum 

-0.25 

2 

2.991 TT 1 0.930 1 2.260 j -0.177 


TABLE I. Locally stable BB soliton pairs for N\ — 1, N 2 ~ 1, 
71 = — 1 , 72 = —1, Vbi = —0.5, Vb 2 = —0.25, belonging to the 
classes BB„ with n = 0,1, 2. The optimal parameters Ax, a\ 
and 02 (units of are obtained by minimizing the global 
potential If (units of 2Er). 


where the Lagrangian density reads 



-Voi cos(2x)|V'iP - ^ cos(2a;)|V'2|'‘ 
-7i2i^in^2p. (17) 


We generalize the Gaussian Ansatz to 


= Aj exp 


. dxoj 
' dt 


(x 


Xoj) + icfijix — Xoj)'^ + idj 


X exp 




(18) 


FIG. 4. Numerical minimization of the pseudopotential 11 [see 
([2^-(23al with soliton peaks around (xoi,xo 2 ) = (0,0) (top 
line), (xoi,xo 2 ) = (tt, 0 ) (central line), (xoi, 2 :o 2 ) = (tt, —tt) 
(bottom line). Plots in the left column represent the varia¬ 
tion of the optimal soliton widths ai (solid blue lines) and 02 
(dashed red lines) with 712 . Plots in the right column show 
the behavior of one optimal width with respect to the other. 
The jumps observed in the central and bottom lines (right 
column) are due to the disappearance of local minima with 
Ax ~ TT and Ax ~ 27r as 712 decreases under a critical value. 
Minima found below the critical 712 coincide with the values 
in the top plots and are thus not shown in the right column. 


(n) (n) 

value 7 cr , such that for 712 < 7 cr metastable solutions 
in BBn would no longer exist. 


In the following, we will numerically analyze the exis¬ 
tence of split BB soliton pairs, as well as their energetic 
and dynamical behavior. To this end, we shall general¬ 
ize our Ansatz to include the positions of the component 
density centers as free parameters. We will also consider 
time-dependent parameters to investigate the dynamics 
of the system. The Gross-Pitaevskii equations Q can be 
restated as a variational problem |49] 


S 


C 



dipj dtp* dtpj dtp*\ 


dxdt = 0 , (16) 


where the variational parameters {Aj,aj,xoj,(pj,6j) for 
j = 1,2 are generally time-dependent and represent, re¬ 
spectively, the amplitude, width, center-of-mass position, 
frequency chirp and overall phase of a soliton in the j-th 
component. The Lagrangian L for the trial wave func¬ 


tions (18), obtained by integrating the Lagrangian den¬ 
sity in (17), reads 


^ + 00 


L = 



— -|- ■ 




- 2a,A= 




\ dt 


_L 4„ . , 2 3 /.2d(('j 

+ 4a,A^- +2a^A^ 


dVo^e ^^cliA‘^ cos(2xoi) 
\/2Vq2G ®"^^a22l2 COs(2xo2) 

_ (a^ni -3:02)^ \ 

47;^2 e “G+< 32 ^ 01^^02342 I 


+ 


-I- 02^ 


(19) 


Note the dependence of the interaction term on Ax, de¬ 
fined in (15). 


The functional derivatives of L with respect to the vari¬ 
ational parameters yield a set of Euler-Lagrange equa¬ 
tions. After appropriate manipulation, we can obtain a 
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FIG. 5. Left panel. Peak positions xoi (solid blue line) and 
(*02) (dotted red line) of optimal BB profiles, obtained by 
searching the local minimum around (aioi, 3 : 02 ) = (tt, —vt), for 
different values of 712. The sudden jump towards (tt, tt) is due 
to the disappearance of the local minimum (see also Fig. 
bottom panels). Right panel. Local minimum value of the 
effective potential energy in Eq. \ 22 \. 


picture for the dynamics of the soliton pairs. In particu¬ 
lar, the equations 


— = 0 , 




dt 


-0 


( 20 ) 

( 21 ) 


can be interpreted as dynamical constraints on the ampli¬ 
tude (particle numbers conservation) and the frequency 
chirp. Taking into account ( 20 )- (ig, the equations of 
motion for the centers of mass and the width can all be 
derived from the effective potential 


n(oi,a2,a;oi,a;o2) = IIi (ai,a;oi) + 112 ( 02 , 0 : 02 ) 
+ ni2 (tti, Xoi, 02,0:02) , 


with 


112 (02,0:02) — 


Ni 

1 N? 

4 al 


+e“ 

cos(2a;oi 

N2 

1 72 iV| 

4 a| 

2 \f^a 2 


-2“2^iV|Vo2 cos(2a;o2) 


ly/ 27 ri 


7 ra 2 


( 22 ) 


( 23 a) 


( 23 b) 


ni2 (oi, Xoi, O2, 0:02) — 


_ (j^oi -X 02 ) 

e “++“2^ 7i2A^iiV2 


( -^12:01 N 
N2X02 

_ 

/^Xoi \ 

^Xo2 I 

(A'i/ 2 )ai 

■■ 

5 a, 

\{N 2 / 2 )a 2 j 


\da2J 


( 23 c) 


+ 02^ 

through the Newton-like equations 


dt^ 


The problem of finding stationary solutions within the 
Ansatz ( 18 ) thus reduces to searching the minima of the 


effective potential. The equations ( 24 ) can also be lin¬ 
earized around the stable equilibrium points for H, to 
obtain information on small center-of-mass and width 
oscillations. On the other hand, due to the (still re¬ 
strictive) form of the trial wave functions, the far-from- 
equilibrium dynamics of Eqs. ( 24 ) cannot be considered 
physically relevant (see the discussion below.) In Table 
[Tj we show some illustrative results of the optimal pa¬ 
rameters Aa:, ai and 02 for solitons in the classes BBn 
with n = 0 , 1 , 2 . The numbers of particles are hxed to 
Ni = N2 = 1 . At 7 i2 = 0 , the three minima are degen¬ 
erate. At 712 = — 0 . 25 , the split solitons in the family 
BB2 are locally stable [see Eq. 0]. while the global 
minimum is, as expected, the overlapped configuration. 
Instead, no local minimum can be found in the family 
BBi, indicating that jcP > — 0 . 25 . 

A more systematic picture of the existence and be¬ 
havior of solitons for small n with varying 712 is given 
in Fig. a numerical minimization procedure is used 
to find stable configurations with their centers of mass 
close to some initial points (xoi,xo2). A minimum with 
centers of mass around (0:01,^02) = (OjO) can be found 
for all negative 712 (top panels). On the other hand, 
the search of minima whose centers of mass are close to 
(xoi,xo2) = (tt, 0) and (ioi,io2) = (-7r,7r) yields dis¬ 
continuous behaviors in the optimal parameters (central 
and bottom panels). The discontinuities are present be¬ 
cause solitons in BBi (respectively BB2) exist only for 

712 > 7cr^ (respectively 712 > 7cr^), while for more neg¬ 
ative values the algorithm actually finds the global min¬ 
imum belonging to BBq. It is also worth noticing that 
in Fig. the optimal widths oi and 02 of the overlapped 
solitons decrease with [712 1 . In the case of split solitons, 
the amplitudes remain almost constant, with a slight in¬ 
crease (more evident in 02) with I712I, due to the attrac¬ 
tion exerted between densities. 

These findings are corroborated by the behavior of the 
center-of-mass positions displayed in Fig. where the 
displacement of the center of mass of the second species 
with increasing interspecies interaction is observed. The 
situation is the same as that depicted in Fig. bottom 
panels. In the left panel of Fig. the jump of 0:02 as 
[7121 is decreased signals the disappearance of the local 
energy minimum. The value of the local minimum of the 
effective potential energy (221 is shown in the right panel. 

In order to check the existence and stability of BB 
soliton pairs as approximate solutions of the GPEs, we 
employ a numerical simulation of the dynamics generated 
by (lU). First, we have checked the stationarity of over¬ 
lapped soliton pairs, localized around a: = 0 . It is possi¬ 
ble to verify, for different values of 712, that the soliton 
pair determined by the minimization procedure is sta¬ 
tionary within very good approximation. In Fig. the 
time evolution of the overlapped solitons is represented 
for 712 = — 0.25 (left) and 712 = —1 (right). Then, we 
have tested the behavior of split soliton pairs in BB2 
in different regimes. In the case 712 = — 0 . 25 , which is 
larger than the critical value 7cr — — 0 . 4 , the split con- 
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FIG. 6. Time evolution of density profiles for overlapped soli- 
tons in the first (blue solid lines) and second (red dashed 
lines) species for 712 = — 0.25 (left panel) and 712 = —1 (right 
panel). The results are obtained by direct numerical integra¬ 
tion. 





FIG. 9 . Split BB solitons inside the regions n = 1 and n = 
2 of the left panel of Fig. at 712 = 0.5 (left panel) and 
712 = 1.0 (right panel). The blue (curve centered at zero) and 
red prohles refer to hrst and second component, respectively, 
while the black line shows the periodicity of the optical lattice. 


IV. SPLITTING DEPENDENCE ON 
INTERSPECIES INTERACTION 


FIG. 7 . Time evolution of density profiles for split solitons 
in the hrst (blue solid lines) and second (red dashed lines) 
species for 712 = — 0.25 (left panel) and 712 = — 0.38 (right 
panel). In the latter case, numerical integration shows the 
instability of the soliton pair. 


figuration evolves in time with slight distortions, but it 
preserves the qualitative features of the initial state for 
all the time of the simulation (left panel of Fig. [^. When 

( 2 ) '—' 

712 — 7cr7 the energetic instability of the soliton pair in 

BB2 is reflected by a dynamical instability: the second- 
species density distribution is gradually attracted by the 
first species (right panel of Fig. [^, ending with an over¬ 
lapped configuration, which is eventually stabilized by 
radiating wave packets [sniisi]. 



FIG. 8. Dependence of chemical potentials (left panel) and 
distance between peaks (right panel) on repulsive interspecies 
interactions of a BB soliton with W = A2 = 1 . Other pa¬ 
rameters are hxed as 71 = 72 = — 1 , Vbi = — 0 . 5 , V02 = — 0 . 25 . 
Top (blue) and bottom (red) curves in the left panel refer 
to chemical potentials of first and second component, respec¬ 
tively, while vertical lines separate 712 regions for existence of 
BB solitons in BBq (overlapped), RRi and BB2 (split). 


In the previous section we observed that split BB soli¬ 
tons can become unstable at some negative critical val¬ 
ues of the interspecies scattering length. We shall now 
investigate these critical values in more detail by direct 
numerical integration of the GPEs, both for attractive 
and repulsive interatomic interactions. 

Let us first discuss the repulsive case. We can con¬ 
sider as initial state an overlapped BB soliton, centered 
at a: = 0, with no interspecies coupling (712 = 0), and 
adiabatically switch on a repulsive interspecies interac¬ 
tion between components at t > 0 . One expects that, 
due to the repulsive interspecies interaction, the initial 
BBq soliton will evolve into a split one belonging to the 
BBi family at some value 712 = 7rep and then into the 


(2) 

BB2 family at 712 = 7rep. This picture coincides with the 
numerical results in Fig.[^ both in terms of the chemical 
potentials and the distances between peaks Ax, normal¬ 
ized to TT. The jumps in the distance are correlated with 
jumps in the chemical potentials at the critical values, 
which are uniquely fixed by the parameters of the system. 
In Fig. 1 ^ the profiles of the split solitons with n = 1 and 
n = 2 are represented, at two different 712 values belong¬ 
ing to their existence curve. Despite the smaller value of 
712, the two BB components in the n = 1 case appear 
to be more distorted in their overlapping region than in 
the n = 2 case. This is an evident consequence of the 
exponential decay of the soliton-soliton interaction with 
distance [see Eq. ( 19 )]. Notice that, as in the attractive 
case, the normalized distance between soliton centers is 
not an integer number. This is a clear consequence of the 
existence of a repulsive force between components. 

When attractive interactions 712 < 0 are adiabatically 
turned on at t > 0, we expect that an initial split soli¬ 
ton in BBriQ, with no > 0, will undergo only one jump 
towards n = 0 . Indeed, from the analysis in the pre¬ 
vious section, we can deduce that the negative critical 
values are ordered as 7^^ > Thus, if 712 > 'yci°\ 





















FIG. 10. Dependence on attractive interspecies interaction of 
chemical potentials (left panel) and distance between density 
centers, xi,X 2 , (right panel) of a BB soliton with A^i = N 2 = 
1. Bottom (blue) and top (red) curves in the left panel refer to 
the first and the second component, respectively, while dotted 
and continuous lines (left panel) and square and triangle sym¬ 
bols (right panel) refer to BB solitons with centers initially 
separated by tt and by 2ti, respectively. Other parameters are 
fixed as 71 = 72 = —1, Vbi = —0.5, V 02 = —0.25. 


the first component is loaded in a linear optical lattice 
and the second component in a nonlinear optical lattice. 
In particular, the existence and stability of split and over¬ 
lapped BB solitons has been investigated by VA, by di¬ 
rect numerical integrations of the coupled GPEs, and by 
direct numerical integrations of the system. The depen¬ 
dence of the existence ranges of BB solitons on the inter¬ 
species interaction parameter has been also investigated. 
In particular, for repulsive interspecies interactions we 
showed the existence of a series of critical values of 712 at 
which transitions from the n- to the u-l-l- split BB soliton 
occur. For attractive interspecies interaction we showed 
that only direct transitions from a split BB solitons to the 
overlapped BB soliton are possible. Since critical values 
at which transitions occur depend on physical parame¬ 
ters of the mixture, these phenomena suggest that split 
BB solitons could be used for indirect measurements of 
these parameters in experiments. 


interactions give enough energy to overcome all the in¬ 
termediate barriers from the UQ-th down to n = 0 . This 
intuitive result, based on energetic considerations, match 
very well the results of the numerical simulation, as one 
can see from Fig. 10 where the cases no = 2 and no = 1 
are represented. 

Since the critical values of 712 at which the transitions 
occur are uniquely fixed by the parameters of the mix¬ 
ture, including the number of atoms and intraspecies in¬ 
teractions, an experimental implementation of the above 
numerical simulations could be used for indirect mea¬ 
surements of the interspecies scattering length of BEG 
mixtures. The interspecies scattering length can also be 
measured from the oscillatory motion of coupled solitons 
as predicted in [ 5 ^ . 


V. CONCLUSIONS 

We have considered matter-wave bright-bright solitons 
in coupled Bose-Einstein condensates, by assuming that 
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